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High- multiplicity final states and transverse-momentum 
dependent parton showering at hadron colliders 
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If large-angle multigluon radiation gives significant contributions to parton showers in 
the LHC high-energy region, appropriate generalizations of parton branching methods are 
required for Monte Carlo simulations of exclusive high-multiplicity final states. We discuss 
the use in this context of transverse-momentum dependent kernels which factorize in the 
region of high energies. We give examples based on ep and pp multi-jet data, and point to 
possible developments for distributions associated with massive final states at the LHC. 
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1 Introduction 

Complex final states with high particle multiplicity are central to many aspects of the LHC 
physics program. Theoretical predictions for these processes require advanced QCD calcu- 
lational tools, which rely both on pcrturbativc results (at present, mostly next-to-lcading- 
order) [T] and on parton shower event generators for realistic collider simulations 

This article discusses aspects of spacelike parton showers that depend on the structure of 
QCD multiparton matrix elements in the multiple-scale region of large center-of-mass energy 
^/s and fixed transferred momenta, and are likely to affect the form of the final states at high 
multiplicity. 

Let us recall that the physical picture underlying the most commonly used branching Monte 
Carlo generators [H [3] is based on coUinear evolution of jets developing, both "forwards" and 
"backwards", from the hard event [3], supplemented (in the case of certain generators) by 
suitable constraints for angularly-ordered phase space [5] . The angular constraints are designed 
to take account of coherence effects from multiple soft-gluon emission [51 [6l [7] . 

The main new effect one observes when trying to push this picture to higher and higher 
energies is that soft-gluon insertion rules [SI [7] based on eikonal emission currents [HI [H] are 
modified in the high-energy, multi-scale region by terms that depend on the total transverse 
momentum transmitted down the initial-state parton decay chain [TUl (TTJ [T^] . As a result, the 
physically relevant distribution to describe initial-state showers becomes the analogue not so 
much of an ordinary parton density but rather of an "unintegrated" parton density, dependent 
on both longitudinal and transverse momenta^ 

The next observation concerns the structure of virtual corrections. Besides Sudakov form- 
factor effects included in standard shower algorithms [U [3] , one needs in general virtual-graph 

^Theoretical aspects of unintegrated pdfs from the point of view of QCD high-energy factorization are 
discussed in 1131 . Associated phenomenological aspects are discussed in 1141 1151 . and references therein (see 
also |16l I17j for recent new work). See works in |18j for first discussions of a more general, nonlocal operator 
formulation of u-pdfs applied to parton showers beyond leading order. 
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terms to be incorporated in transverse-momentum dependent (but universal) splitting func- 
tions [10l[T8l[T9l[20l[2T] in order to take account of gluon coherence not only for coUinear-ordered 
emissions but also in the non-ordered region that opens up at high ^fsj-p^. 

These finite-k_L corrections to parton branching have important implications for multiplicity 
distributions and the structure of angular correlations in final states with high multiplicity. In 
the next section we discuss examples of such effects in the case of multi-jet production in ep 
and "pp collisions. In Sec. 3 we go on to possible developments involving the hadroproduction of 
massive states. In particular, we point to studies beginning to investigate the role of showering 
corrections versus multiparton interaction corrections in Monte Carlo event generators [551 HSl 
We give final remarks in Sec. 4. 

2 Jet-jet correlations 

In a multi-jet event the correlation in the azimuthal angle A(/) between the two hardest jets pro- 
vides a useful measurement, sensitive to how well QCD multiple-radiation effects are described. 
In leading order one expects two back-to-back jets; higher-order radiative contributions cause 
the A(/) distribution to spread out. Near A(/) ^ tt the measurement is mostly sensitive to in- 
frared effects from soft-gluon emission; the behavior as Acjj decreases is driven by hard parton 
radiation. At the LHC such measurements may become accessible relatively early and be used 
to test the description of complex hadronic final states by Monte Carlo generators [25] . 

Experimental data on Acf) correlations are available from the Tevatron [26j and from Hera [27l 
[28| . These analyses indicate that the comparison of data with Monte Carlos and perturbative 
results are very different in the two cases. Observe in Fig. [T] that the Tevatron A(/) distribution 
drops by about two orders of magnitude over a fairly narrow range, essentially still close to 
the two-jet region. The measurement is dominated by leading-order processes, with small 
sub-leading corrections. Correspondingly, data are reasonably well described both by coUinear 
showers (Herwig and new Pythia tuning) and by fixed-order NLO calculations [25], [26]. 

The Hera A^ measurements, on the other hand, are much more sensitive to higher orders. 
Fig. [2] [27]. NLO results for di-jet azimuthal distributions are affected by large corrections in 
the small-A(/> and small- a; region, and begin to fall below the data for three-jet distributions in 
the smallest A<j> bins . These measurements are likely relevant for extrapolation of initial- 
state showering effects to the LHC, given the large phase space available for jet production, 
and relatively small ratio of jet transverse energy to the center-of-mass energy. 

Refs. |29l [30] analyze the effects of finite-kj^ corrections to initial-state showers, using 
data |27j on jet angular and momentum correlations, and factorization at fixed transverse 
momentum [11] valid for high energies. Fig. [3] shows results from the coUinear Herwig Monte 
Carlo [3T] and from the kj^-shower Cascade Monte Carlo [K] for the distributions in A(f> and 
Apt |27l 129] . measuring the transverse momentum imbalance between the leading jets. The 
largest differences between the two Monte Carlos arc at small A(f> and small Apt, where the 
two highest Et jets arc away from the back to back region and one has effectively three hard, 
well-separated jets. By examining the angular distribution of the third jet, Ref. [29] finds sig- 
nificant contributions from regions where the transverse momenta in the initial state shower 
are not ordered. The description of the measurement by the kj^-shower is good, whereas the 
collinear-based Herwig shower is not sufficient to describe the observed shape. 

The physical picture underlying the k^-shower method involves both transverse-momentum 
dependent pdfs and matrix elements. Fig. [4] [30] illustrates the relative contribution of these 
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Figure 1: Dijet azimuthal correlations measured by DO along with the Herwig and Pythia 
results pS] , 



different components to the result, showing different approximations to the azimuthal dijet dis- 
tribution normalized to the back-to-back cross section. The solid red curve is the full result 
The dashed blue curve is obtained from the same unintegrated pdf 's but by taking the coUinear 
approximation in the hard matrix element. The dashed curve drops much faster than the full 
result as decreases, indicating that the high-k_L component in the ME [5D] is necessary to 
describe jet correlations for small A0. The dotted (violet) curve is the result obtained from the 
unintegrated pdf without any resolved branching. This represents the contribution of the in- 
trinsic distribution only, corresponding to nonperturbative, predominantly low-k_L modes. That 
is, in the dotted (violet) curve one retains an intrinsic k^ ^ but no effects of coherence. We 
see that the resulting jet correlations in this case are down by an order of magnitude. The 
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Figure 2: (left) Azimuth dependence and (right) Bjorken-x dependence of ep di-jet distribu- 
tions [23, compared with NLO results. 
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Figure 3: (left) Angular correlations and (right) momentum correlations [29] in three-jet final 
states measured by (27] j compared with kj^-shower (Cascade) and coUincar-shower (Herwig) 
Monte Carlo results. 



inclusion of the pcrturbatively computed high-kx correction distinguishes the calculation |29| 
from other shower approaches (see e.g. [32]) that include transverse momentum dependence in 
the pdfs but not in the matrix elements. 
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Figure 4: The dijet azimuthal distribution [30j normalized to the back-to-back cross section: 
(solid red) full result (u-pdf ME); (dashed blue) no finite-k^ correction in ME (u-pdf 
MEco/hri.); (dotted violet) u-pdf with no resolved branching. 

The above observations underline the role of accurate multi-jet measurements in events 
associated with proton scattering off virtual photons [16j . Further phenomenological analyses 
of available jet correlation data (and multiplicity distributions [29]) will be helpful. In this 
respect note that the effects of coherence emphasized above dominate for sufficiently small 
and small x. But the unintegrated formulation of parton showers is potentially more 
general [TH [331 [31] . Such analyses can be of use in attempts to relate [35] shower effects in DIS 
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event shapes [36] measuring the transverse momentum in the current region to vector-boson 
hadroproduction pr spectra [22\ [37] . 

3 Massive final states 

Corrections to collinear-ordered showers affect heavy mass production, including the structure 
of the final states associated with heavy flavor and heavy boson production. We next point to 
examples that depend on the physics of unintegrated gluon distributions. 

Measurements of angular correlations for bottom quark jets have recently been performed at 
the Tevatron [351 [SSI iO]- See [T1[1T1|11[13] for reviews of related phenomenology Results for b- 
jet distributions in invariant mass and azimuthal angle are shown in Fig. [51 [55] and Fig. [10] . 
Monte Carlo simulations based on Pythia, Herwig and MC@NLO do not appear to give 
satisfactory descriptions of the observations [U [32] especially at small Acj). The measurement of 
&-jet correlations has considerable interest, given their potential sensitivity to soft underlying 
events [39l|44] and possibly models for multiple-parton interactions [221[23]lj In this context it is 
worth noting the possible role of showering corrections, at the level of single-parton interactions, 
due to transverse-momentum dependent parton branching. 
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Figure 5: Invariant-mass distribution and azimuthal- angle distribution for production of &-jcts 
at the Tevatron [39] . 

To this end let us recall that heavy flavor hadroproduction is dominated for sufficiently 
high energies by gluon splitting into heavy-quark pairs [T3j, g — > QQ where g is produced 
from the spacelike jet. The high-energy asymptotic behavior is controlled by a triple-pole 
singularity |13j in the complex plane of the Mellin moment conjugate to the transferred k_L. 
The coefficient functions associated to this singularity enhance regions that are not ordered in 
transverse momentum in the initial state shower. In fact, such contributions are already found 
to be significant at the level of the NLO correction [42l [43] H 

^This is unlike the DIS jet correlation data discussed in the previous section, where multiparton contributions 
are believed to be much suppressed |16) . 

^It is possible that terms of this kind at orders higher than NLO are responsible for the rather large theoretical 
uncertainties found |42l I43j in the NLO predictions when going from the Tevatron to the LHC. 
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Figure 6: Comparison of data and Monte Carlos for &-jet azimuthal correlations at the Teva- 
tron [40] , 

A typical contribution to g QQ is pictured in Fig. [7^. Note that for small this graph 
gives effectively a contribution of leading order. Corrections of the next order from additional 
jet emission are shown in Fig.[7j5. In the notation of Fig. [71 the triple-pole behavior is produced 
from regions in which TOg <C (fcr + ^t)^ ^ fc|i ~ fc^, where tjiq is the heavy quark mass, and 
Rt, kip are transverse momentum vectors. Collinear shower calculations, even if supplemented 
by NLO matrix elements as in MC@NLO, are not designed to take account of this behavior. 
This is likely to reduce the numerical stability of predictions as one goes to higher and higher 
energies. It may cause a non-negligible contribution from showering to be missed in the 6-jet 
A(j) distribution at small A(f>. On the other hand, such corrections can be obtained by methods 
based on transverse-momentum dependent parton branching, as those discussed in the previous 
section. It is of interest to analyze these contributions in comparison with those, e.g. in Fig. [SJ 
from multiple interactions. 



Figure 7: (a) Heavy quark hadroproduction from gluon showering at high energy; (b) next 
correction from extra jet emission. 

Such an analysis would pay off as one goes from the Tevatron to the LHC. Parton showers 
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at the LHC will be more influenced by the asymptotic high-energy pole. Let us observe that 
effects of a similar physical origin will affect the structure of final states associated to production 
processes predominantly coupled to gluons, e.g. central scalar boson production [45]. See studies 
of showering effects in this case [53]. These can affect the description of soft underlying events 
and minijets [44l[22] as well as the use of exclusive production channels [46] . 

4 Concluding remarks 

Final states with high particle multiplicity acquire qualitatively new features at the LHC com- 
pared to previous collider experiments due to the large phase space opening up for events 
characterized by multiple hard scales, possibly widely disparate from each other. This brings 
in potentially large perturbative corrections to the hard-scattering event and potentially new 
effects in the parton-shower components of the process. 

If large-angle multigluon radiation gives significant contributions to parton showers at the 
LHC, appropriate generalizations of parton branching methods are required. We have discussed 
applications of transverse-momentum dependent kernels for parton showering that follow from 
factorization properties of QCD multiparton matrix elements in the high energy region, which 
are valid not only for coUinear emission but also at finite angles. 

While we have focused on observables that are sensitive primarily to the physics of initial- 
state gluonic showers, expressible in terms of "unintegrated" gluon densities, treatments of 
quark contributions to showers at unintegrated level are being worked on (see e.g. |16l I17|V 
In this respect, theoretical results for splitting kernels [20] already applied to inclusive phe- 
nomenology could also be of use in calculations for exclusive final states. 

We have considered jet and heavy-mass production processes in the central rapidity region. 
Techniques are being developed [17] to allow one to also address multi-particle hard processes 
at forward rapidities. 
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